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ABSTRACT 

Some recent experiments detecting very high energy (VHE) gamma-rays above 10-20 TeV indepen- 
dently reported VHE bursts for some of bright gamma-ray bursts (GRBs). If these signals are truly from 
GRBs, these GRBs must emit a much larger amount of energy as VHE gamma-rays than in the ordinary 
photon energy range of GRBs (keV-MeV). We show that such extreme phenomena can be reasonably 
explained by synchrotron radiation of protons accelerated to ~ io 20 ~ 21 eV, which has been predicted 
by Totani (1998a). Protons seem to carry about (m p /m e ) times larger energy than electrons, and hence 
the total energy liberated by one GRB becomes as large as ~ 10 56 (Af2/47r) ergs. Therefore a strong 
beaming of GRB emission is highly likely. Extension of the VHE spectrum beyond 20 TeV gives a nearly 
model-independent lower limit of the Lorentz factor of GRBs, as 7 > 500. Furthermore, our model gives 
the correct energy range and time variability of ordinary keV-MeV gamma-rays of GRBs by synchrotron 
radiation of electrons. Therefore the VHE bursts of GRBs strongly support the hypothesis that ultra 
high energy cosmic rays observed on the Earth are produced by GRBs. 

Subject headings: acceleration of particles — cosmic rays — gamma rays: bursts — gamma rays: theory 



Two major ground-based cosmic ray/gamma-ray detec- 
tors, the Tibet air shower array (Amenomori et al. 1996) 
and the HEGRA AIROBICC Cherenkov array (Padilla 
et al. 1998) have independently reported significant ex- 
cesses of 10-20 TeV gamma-rays coincident with some 
GRBs both in direction and burst time. There are 57 
GRBs detected by the BATSE (Meegan et al. 1996) in 
the field of view of the Tibet array, and the Tibet group 
found that several GRBs out of the 57 GRBs show sig- 
nificant excess of 10 TeV gamma-ray events with a time 
scale of ~ 10 seconds. The statistical significance was es- 
timated to be about 6 sigma. There are four GRBs in 
the field of view of the HEGRA array, and the HEGRA 
group observed 11 gamma-ray like events above 20 TeV 
in a 4-minutes bin coincident with GRB920925c in the 
GRANAT/WATCH GRB catalog (Sazonov et al. 1998), 
whereas only 0.93 events are expected as background (5.4 
sigma). The chance probability taking into account an ap- 
propriate trial factor was estimated as 0.3 %. Seven out 
of the 11 gamma-ray events are clustered within 22 sec- 
onds, suggesting that the burst time scale is ~ 10 seconds 
also for this GRB. Considering the fact that two different 
experimental groups independently reported similar sig- 
nificant signals of 10-20 TeV gamma-rays, now we must 
seriously take into consideration these very high energy 
(VHE) emissions. Although further observations are nec- 
essary to confirm these indications of VHE bursts, in this 
letter we show that these VHE bursts, if confirmed, would 
result in a drastic change of our picture of GRBs. 

There is a serious problem in interpreting these VHE 
gamma-rays. The cosmological origin of GRBs has al- 
most been confirmed (Metzger et al. 1997; Kulkarni et 



al. 1998) and most of GRBs are considered to come from 
large distances beyond redshift z ~ 1. It is well known 
that VHE gamma-rays above ~ TeV from such cosmologi- 
cal distances are not visible because they suffer serious at- 
tenuation by electron-positron pair creation with the inter- 
galactic infrared background radiation (Stecker, de Jager, 
& Salamon 1992 and references therein). A recent cal- 
culation (Stecker & de Jager 1998) shows that the optical 
depth r for 10 TeV gamma-rays is about 5-10 (correspond- 
ing to an attenuation factor of e~ r ~ 7 x 10~ 3 -5 x 10 -5 ) 
for the source redshift of z = 0.1-0.3, and hence there 
is a horizon for such gamma-rays at z ~ 0.2. Since the 
GRB920925c has a large energy fluence (Sazonov et al. 
1998) of 2.1 x 10~ 5 erg cm" 2 in 20-60 keV (typical fluence 
range of the BATSE catalog is 10- 7 -10~ 4 erg cm 2 in 50- 
300 keV), it is likely that the GRBs having the VHE emis- 
sion are nearby bursts within z < 0.2. In fact, with typical 
distance scales of z max = 1-3, where z max is the redshift 
of the most distant GRBs observed by the BATSE, a few 
percent of the BATSE GRBs are expected to be located 
around z ~ 0.1-0.3. This is consistent with the fact that 
only a few out of 57 BATSE GRBs in the Tibet data seem 
to have VHE emission. It should be noted that, although 
the GRBs having the VHE emission must be the nearest 
bursts, the flux of keV-MeV gamma-rays is not necessar- 
ily the brightest because there may be intrinsic luminosity 
dispersion of GRBs. 

The energy fluence of 10 TeV gamma-rays from 
GRB920925c was estimated (Padilla et al. 1998) as 
6.9 x 10 -6 erg cm" 2 , which is comparable to that in the soft 
gamma-ray range of 20-60 keV (2.1 x 10 -5 erg cm -2 ). The 
VHE fluence of a few GRBs observed by the Tibet array 
was estimated as 10 _6 -4 x 10 -5 erg cm -2 (Amenomori et 
al. 1996), which is also comparable to that of soft gamma- 
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rays of relatively bright GRBs. Furthermore, we must 
take into account the fact that the intergalactic attenu- 
ation factor of 10-20 TeV gamma-rays is likely as small as 
- lO^-LO" 3 for z ~ 0.1-0.2. This suggests that the en- 
ergy emitted as VHE gamma-rays is much larger than that 
emitted as keV-MeV gamma-rays by a factor of ~ 10 2 - 
10 3 . This emission is extraordinarily energetic, but in the 
following we show that such VHE emission can be reason- 
ably explained by synchrotron radiation of ~ 10 20 ~ 21 eV 
protons in GRBs. The origin of ultra high energy cosmic 
rays (UHECRs), whose energy extends beyond 10 20 eV, 
is one of the most interesting mysteries in astrophysics as 
well as GRBs, and the observations of 10 TeV gamma- 
rays by the Tibet and HEGRA groups strongly suggest 
that these two mysterious phenomena have a common ori- 
gin, as it has been speculated (Waxman 1995; Milgrom & 
Usov 1995; Vietri 1995). 

Since the energy emitted as 10 TeV gamma-rays is much 
larger than that of soft gamma-rays, the VHE emission is 
likely produced by protons while the keV-MeV emission is 
produced by electrons. The current theory of GRBs (e.g., 
Piran 1997) attributes the origin of gamma-ray emission to 
dissipation of the kinetic energy of ultra relativistic bulk 
motion with a Lorentz factor of 7 ~ 10 2 -10 3 , in inter- 
nal shocks generated by relative difference in speeds of 
ejected matter. Therefore at least in the initial stage pro- 
tons must carry much larger energy than electrons by a 
factor of ~ m p /m e ~ 2, 000. The emission mechanism of 
keV-MeV gamma-rays is generally considered to be elec- 
tron synchrotron. If the energy carried by protons is ef- 
fectively emitted as VHE photons before the equipartition 
between protons and electrons is achieved, the very strong 
flux of 10 TeV gamma-rays observed by the two detectors 
is naturally explained. Since the energy emitted as keV- 
MeV gamma-rays is typically 10 52-53 erg for cosmological 
GRBs, the total energy liberated by one GRB event must 
be as large as ~ 10 56 erg if radiation is isotropic. This en- 
ergy may seem extremely large, but if GRBs are strongly 
beamed, such a huge amount of energy can be supplied. 
For example, the 'microquasar model' (Paczynski 1998) 
can produce about 5 x 10 54 erg and the above energy can 
be explained if the beaming factor, 47r/Af2, is larger than 
~ 20. Such beaming factor is not unreasonable because 
the origin of energy of the microquasar model is rotation 
energy of a Kerr black hole with a mass of ~ 10 M Q . Ra- 
diation likely deviates from spherical symmetry because of 
the rotation. 

In fact, we have pointed out (Totani 1998a) that the en- 
ergy transfer between protons and electrons may be a quite 
inefficient process and strong afterglow of GRBs in TeV 
range may exist by synchrotron radiation of protons accel- 
erated up to 10 20 " 21 eV. The VHE bursts beyond 10 TeV 
suggested by the Tibet and HEGRA experiments may also 
be explained by the proton synchrotron radiation. How- 
ever, in the previous work we have considered the after- 
glow phase, i.e., the external shock generated by collisions 
of ejected matter and ambient interstellar matter. In this 
case, the cooling time scale of such protons is about a 
few days and hence the afterglow model cannot explain 
the observations of the two detectors showing a burst-like 
time structure of ~ 10 seconds. These VHE emissions are 
coincident in time with keV-MeV GRBs (Amenomori et 



al. 1996; Padilla et al. 1998), and should be explained 
by internal shocks rather than external shocks. Here we 
present an internal shock model which can explain the ob- 
served VHE emission by proton synchrotron as well as the 
ordinary keV-MeV emission by electron synchrotron. 

Suppose that the total energy E is isotropically emitted 
during a time interval T = lOTio sec which is the typi- 
cal duration of GRBs, and this energy is converted into 
kinetic energy of relativistic bulk motion with a Lorentz 
factor of 7 = 3007300- As discussed above, we consider a 
typical case with E — 10 56 i?56 erg. The following analysis 
is not affected by the unknown beaming factor if we scale 
the total energy appropriately. If the Lorentz factor has 
fluctuations of A7 ~ 7, faster ejecta will catch up with 
slower ejecta at 

R ~ j 2 T = 2.7 x 10 16 7 2 00 T 10 cm . (1) 

(We use the natural units of c = h — 1.) Therefore 
the energy density of internal shocks can be estimated as 
Piab ~ E/ (AttR 2 T) at the laboratory frame. The restframe 
energy density p rcst ~ puh/l 2 is then given as 

/>rest ~ AlT R2 Tl 2 = 4 X l ^ ^1^^ Crg Cm" 3 . (2) 

We assume that there is a magnetic field which is nearly 
in equipartition with this energy density, i.e., 

B = 3.2 x 10 3 ^ /2 ^56 /2 730 3 o7 i ro 3/2 Gauss , (3) 

where £b is a parameter of order unity for degree of 
equipartition between magnetic fields and shock heated 
matter. 

The most efficient energy loss process of protons or elec- 
trons would be synchrotron radiation. The photon en- 
ergy of proton synchrotron in observer's frame is given 
by e 7 = £ 2 e£?/(m 3 7), where e p is the proton energy in 
the observer's frame. Therefore, if the observed 10 TeV 
gamma-rays are proton synchrotron, the energy of parent 
protons must be 

£p = 3 ' 6 X 1Q2 ° (i^y^W 7 ^ 174 ^ 74 eV , 

(4) 

which is just the energy scale of the UHECRs observed 
on the Earth. It has been known that physical quan- 
tities of GRBs allow protons to be accelerated to such 
high energy when E <~ 10 51 erg (Waxman 1995; Vietri 
1995). We are now considering a much more energetic 
case and we have to check whether protons are acceler- 
ated to ~ 10 20 eV in the present model. The acceleration 
time of Fermi acceleration is roughly given by the gyrope- 
riod, 2nr]rL, where n, — m p ^ p /(eB) is the Larmor radius 
of protons, -f p — e p j (m p j) the Lorentz factor of protons in 
the shock restframe, and r\ a parameter of order unity (e.g., 
de Jager et al. 1996). The acceleration limit is given by 
2-KrirL = R/~f, where R/j is the expansion time measured 
in the shock restframe. (The shell thickness measured in 
the shell frame is also given by ~ T7 ~ R/j, and the con- 
finement condition is also satisfied automatically.) Hence, 
the maximum proton energy is given by 

e p = 4.1 x 10 21 7 7 - 1 7 3 oUB 72 ^5 1 6 72 rr o 172 eV . (5) 
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Proton acceleration is also limited by synchrotron cooling, 
i.e., 2irr]rL = t coo i, where < coo i = 6nmp (a T mlB 2 "/ p ) is the 
synchrotron cooling time at the shock rest. This suggests 
that the maximum energy achieved by proton acceleration 
is 

e p = 4.4 x lO^/^E-^^TU 4 eV . (6) 

This energy is higher than that corresponding to 10 TeV 
synchrotron photons (see eq. |^), and hence proton syn- 
chrotron radiation extends to ~ 10 TeV range. In fact, 
the maximum photon energy of synchrotron emission is 
determined only by fundamental constants independently 
of B, if the maximum particle energy is constrained by 
synchrotron cooling. For electron synchrotron, the max- 
imum photon energy is given by ~ 25?7~ 1 7(sin0)~ 1 MeV 
(de Jager et al. 1996), and that for proton synchrotron 
is higher than this by a factor of (m p /m e ), where 9 is 
the pitch angle of particles along the direction of mag- 
netic fields. By using (sin 6) — ir/4 = 0.785 for isotropic 
particle distribution, the maximum photon energy of pro- 
ton synchrotron is ~ 18?7 -1 73oo TeV. Since this maximum 
energy is quite model-independent and rj is likely larger 
than 1 (de Jager et al. 1996), we can infer the Lorentz 
factor of GRBs from the observed maximum energy of 
VHE gamma-rays. The energy threshold (50% trigger ef- 
ficiency) of the HEGRA AIROBICC array is 16 (25) TeV 
(Padilla et al. 1998), and hence we obtain a lower limit of 
7 > 3007? and likely > 500r?. 

It should be noted that the maximum proton energy is 
determined by the synchrotron cooling constraint rather 
than the expansion time. This suggests that protons 
around the maximum energy will effectively lose their en- 
ergy in GRBs and hence significant fraction of energy car- 
ried by protons will be converted into 10 TeV photons. 
Such VHE emission would be much more energetic than 
keV-MeV emission by electrons, giving a natural explana- 
tion for the VHE bursts observed by the two experiments. 
In fact, the cooling time of protons for an observer is given 

by 

tcool.obs- ^ - U ' D ^ 10TeV J 

X73 4 ooG 3/4 ^ 6 3/4 < 4 sec , (7) 

which is sufficiently short to explain the time structure of 
the observed TeV bursts with the time scale of ~ 10 sec. 

Next we consider whether electrons in this model can 
explain the ordinary keV-MeV GRBs. The photon energy 
of electron synchrotron in observer's frame is given by 

e 7 = ^ = 1.0 fc) 2 C^E^T-r keV , (8) 

where j e is the electron Lorentz factor measured in the 
shock frame. If the energy transfer from protons to elec- 
trons is inefficient, which is the case we are now consid- 
ering, the minimum Lorentz factor of electrons is a direct 
consequence of initial bulk motion and hence 7 e ~ 7. Ob- 
servationally the lower bound of energy range of GRBs is a 
few keV, and our model naturally gives this lower bound. 
Electrons are also accelerated to higher energies and such 



electrons explain the higher energy gamma-rays extending 
up to MeV-GeV. As is done for protons above, we can 
estimate the maximum energy of electrons assuming that 
the acceleration time is 2-Kr]r^. The maximum electron en- 
ergy is determined by the synchrotron cooling constraint 
and hence the maximum photon energy of synchrotron is 
given by the formula: 

£-« = 25 rf 1 1 (sin^r 1 MeV (9) 
= 9.5 ttSoo GeV . (10) 

This is consistent with the observation of some bright 
GRBs extending to the GeV range (Hurley 1996). [Pro- 
ton synchrotron may also contribute (Totani 1998a) to the 
GeV range gamma-ray emission, especially for the long du- 
ration GeV emission observed for GRB940217 (Hurley et 
al. 1994).] The cooling time of electrons for an observer 
is written as: 

Wob , = ^ = u 6x io-.y^)- 1/2 

^W'-^t'X' 4 msec , (11) 

which is near the shortest time variability observed in 
GRBs (~ msec). Hence we can interpret the msec time 
scale of variability as synchrotron cooling time of electrons. 
This hypothesis predicts that the time scale of variability 

— 1/2 

changes with photon energy as oc e 7 , which can be 
tested by existing data. 

We have shown that the recent observations of 10 TeV 
gamma-rays from some GRBs can be explained quite well 
if protons carry much larger energy than electrons and 
protons are accelerated up to ~ 10 20-21 eV. Furthermore, 
electron synchrotron radiation in our model naturally gives 
the correct energy range and time variability of the ordi- 
nary keV-MeV gamma-rays of GRBs. It seems quite diffi- 
cult to explain such strong emission of 10 TeV gamma-rays 
by other radiation processes and hence the observations 
of the two independent experimental groups suggest that 
protons in GRBs are actually accelerated up to the en- 
ergy scale of the UHECRs observed on the Earth. The 
flux of UHECRs is difficult to estimate theoretically be- 
cause the escape fraction of such protons from GRBs is 
poorly known and the lifetime of UHECRs in the inter- 
galactic space drastically changes around 10 20 eV. This 
is known as the Greisen-Zatsepin-Kuzmin cutoff (Greisen 
1966; Zatsepin & Kuzmin 1966), which is due to the photo- 
pion production by the interaction between UHECRs and 
the 2.7K cosmic microwave background radiation. The 
lifetime is about tcr ~ 10 8 yrs for protons with energy 
of ~ 2 x 10 20 eV (Berezinskh et al. 1990), corresponding 
to a horizon for UHECRs at ~ 30 Mpc. However, be- 
cause of the large total energy of ~ 10 56 6 _1 erg for one 
GRB, where b = 47r/AJ7 is the beaming factor, it can be 
concluded that the production rate of UHECRs in GRBs 
is sufficiently large to explain the observed UHECR flux 
above 10 20 eV. In fact, we can calculate the energy pro- 
duction of UHECRs per one GRB required to explain the 
observed UHECR flux as 

E CK = 2 x 10^ ( 

\7 x lO^ 1 erg cm -13 / 
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/ t cr \ 1 / Rgrb \ 1 , 19 % 

X ^ 10 8 yr) \ 10-% yr-iMpc-3 J CTg ' [ ' 

where pen is the energy density of UHECRs deduced 
from the observed UHECR flux (Yoshida and Dai 1998), 
and i?GRB the present-day rate density of GRBs (Totani 
1997; 1998b). If the total energy liberated by one GRB is 
- lO 56 ^ 1 erg, UHECR production with E C r ~ lO 54 ^ 1 
erg would be quite reasonable and hence we conclude 



that the 10 TeV gamma-rays observed by the Tibet and 
HEGRA groups give a strong support for the hypothesis 
that UHECRs are produced by GRBs. 
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